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V391 Peg (HS 2201+2610) is an extreme horizontal branch subdwarf B (sdB) star, it is an hybrid pulsator 
showing p- and p-mode oscillations, and hosts a 3.2/sini Mjup planet at an orbital distance of about 1.7 
AU. In order to improve the characterization of the star, we measured the pulsation amplitudes in the u'g'r' 
SLOAN photometric bands using ULTRACAM at the William Herschel 4.2 m telescope and we compared 
them with theoretical values. The preliminary results presented in this article conclusively show that the 
two main pulsation periods at 349.5 and 354.1 s are a radial and a dipole mode respectively. This is the 
first time that the degree index of multiple modes has been uniquely identified for an sdB star as faint as 
V391 Peg (B=14.4), proving that multicolor photometry is definitely an efficient technique to constrain mode 
identification, provided that the data have a high enough quality. 
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1 Introduction 

About half of field sdB stars, which reside in binary 
systems, can form through common envelope ejection 
or stable Roche lobe overflow (Han et al. 2002, 2003). 
It is more difficult to explain the formation of a single 
sdB star. Two scenarios are possible: the merger of two 
low-mass helium white dwarfs and an early hot helium 
flash; but both are not fully consistent with the obser- 
vations. A recent review on these arguments is given by 
Heber (2009). Another possibility, suggested by Soker 
(1998), is that the huge mass loss needed to form an sdB 
star is triggered by low-mass bodies, planets or brown 
dwarfs (BDs). Although this possibility has not been 
tested by detailed models yet, a planet to the pulsat- 
ing sdB star V391 Peg (Silvotti et al. 2007) and three 
circumbinary planets/BDs to the eclipsing sdB-f M sys- 
tems HW Vir (Lee et al. 2009) and HS 0705-1-6700 
(Qian et al. 2009) suggest that sdB planets/BDs could 
be a relatively common phenomenon (see also Bear & 



Soker 2010). A systematic search for sdB substellar ob- 
jects around 4 sdB stars using the timing method is 
the main goal of the EXOTIME project (Schuh et al. 
2010, Benatti et al. 2010). 

V391 Peg (HS 2201-1-2610) is a particularly inter- 
esting system formed by an sdB star and a 3.2/sim 
Mjtip planet orbiting the host star in 3.2 years at a dis- 
tance of about 1.7 AU (Silvotti et al. 2007). The sdB 
star is a hybrid pulsator showing p and g-mode oscilla- 
tions at the same time (0stensen et al. 2001, Lutz et al. 
2009), offering a unique opportunity to characterize the 
host star through asteroseismic methods. A preliminary 
mode identification of the higher pulsation frequencies 
(p-modes) was proposed in Silvotti et al. (2002): the 
two main pulsation periods of 349.5 and 354.1 s could 
be reproduced with /=0 (fc=l) and 1=1 (fc=l) respec- 
tively. However this solution was not unique due to the 
small number of detected frequencies and other solu- 
tions could not be excluded. 



Based on observations obtained at the 4.2 m WHT 
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Fig. 1 Upper panels: representative u'g'r' light curves (20 October 2007), with beating effects clearly visible. 
Lower panels: amplitude spectra of the whole 8-nights run showing the two regions of excited p- and g-modes. 



2 Observations, reduction and analysis 

V391 Peg was observed for 8 consecutive nights, from 
October 16 to 23, 2007, using ULTRACAM (Dhillon 
et al. 2007) at the 4.2 m William Herschel telescope 
(WHT). In total we collected 260,592 exposures in 
three photometric bands {u'g'r') of the SLOAN sys- 
tem, with exposure times between 1.2 and 1.6 s and a 
dead time of only 25 ms between one frame and the 
next. The reduction was performed using the ULTRA- 
CAM pipeline (see, for example, Littlefair et al. 2008), 
including bias and flat field correction and aperture 
photometry. Then we computed differential photom- 
etry, dividing the target's counts by the counts of a 
comparison star, we binned the data to an effective ex- 
posure time of 10 s and we performed a correction for 
the residual extinction. The last step is crucial for the 
g-modes which have periods of '^0.5-2 h and are par- 
ticularly disturbed by variations that occur on similar 
time scales. Finally we applied the bary centric correc- 
tion to the times. More details regarding observations 
and data reduction will be given in a forthcoming pa- 
per (Silvotti et al. in preparation). At the end of the 



process, for each filter we obtained a single file with 2 
columns: barycentric Julian date and fractional bright- 
ness variation. 

These files were analyzed using Period04 (Lenz and 
Breger 2005) in order to determine the pulsation fre- 
quencies and the amplitudes in the different bands. A 
portion of the light curves and the amplitude spectra 
in the three photometric bands are shown in Fig. 1. 
The frequencies obtained were compared with those ob- 
tained from previous runs and indeed we found a per- 
fect agreement except for f4 and f5: the value 2921.816 
/iHz found previously for f4 (Silvotti et al. 2002) is 
now 2910.272 ^Hz, indicating that the old value was 
a 1-day alias of the correct frequency. The new value 
of f4 is confirmed by two other independent runs with 
high frequency resolution at the 3.6 m TNG (August 
2008) and 1.3 m MDM (October 2007) telescopes (Sil- 
votti et al. in preparation). F5 (2738.015 /zHz, Silvotti 
et al. 2002) is not found in any of the new observa- 
tions and the WHT/ULTRACAM data suggest 2 new 
low-amplitude frequencies. An updated list of frequen- 
cies, including also the low-frequency g-modes, will be 
published (Silvotti et al. in preparation). 
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Fig. 2 Fit to the two dominant modes of V391 Peg, including non-adiabatic effects. 



Using the improved list of frequencies, we measured 
the amplitudes of the various frequencies by means of 
least-square sinusoidal fits. In this paper we concen- 
trate only on fl and f2, for which the errors on the 
amplitudes are sufficiently small to obtain significant 
results. 



The monochromatic atmospheric quantities were 
then convolved over the ULTRACAM filters, taking 
into account the filter response curves as well as the 
quantum efficency of the CCDs, and the transparency 
curve for a representative observatory site (we used 
Cerro Tololo which is at a similar altitude as the La 
Palma observatory). 



3 Comparison with theoretical 
amphtudes and mode identification 

The amplitudes obtained have been compared with the- 
oretical amplitudes calculated following the same pro- 
cedure as described in Randall et al. (2005). As in- 
put we took the atmospheric parameters from 0stensen 
et al. 2001 (Toff =29,300 K, log.g=5.4). With the very 
small error bars of the ULTRACAM amplitudes, the 
quality of the fit is very sensitive to the exact input 
values of the atmospheric parameters. However, taking 
into account the T^s logg uncertainties, we never ob- 
tain any change in the mode identification. The same 
is true when we use slitghtly different atmospheric pa- 
rameters obtained by one of us (G.F.) on the basis 
of spectra acquired at the 90 inch telescope at Kitt 
Peak (Teff =30,000 K, logg=5.46) or at the MMT 
(reff=29,900 K, log (7=5.51), made available to us by 
Betsy Green (private communication). 



Non-adiabatic effects were computed as they sig- 
nificantly infiuence the theoretical colour-amplitudes. 
They were estimated using the equations given in Ran- 
dall et al. (2005) from the adiabatic and non-adiabatic 
eigenfunctions calculated from second generation enve- 
lope models (Charpinet et al. 1997), using the follow- 
ing stellar parameters: Toff =29, 300 K, log (7=5.4, total 
stellar mass M*=0.48 Mq and fractional thickness of 
hydrogen-rich envelope log (7(7f)=— 2.95, although the 
two last values do not really influence the results much, 
as long as they take on reasonable values. The model 
was then analysed with adiabatic and non-adiabatic 
pulsation codes, and the non-adiabatic quantities R 
and vPt (defined in Randall et al. 2005) were computed 
for the period spectrum of the model. Since R and ^ t 
are not dependent on the degree index while they de- 
pend quite strongly on the period of the mode in ques- 
tion, their value were interpolated in period space to 
find their optimal values for the observed periods. 
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The wavelength-integrated atmospheric quantities 
and the non-adiabatic parameters were then used to 
calculate the theoretical amplitudes. As a last step, the 
theoretical amplitudes were fit to those observed using 
the x^ minimisation technique described in Randall et 
al. (2005). 

The results, shown in Fig. 2, indicate a unique solu- 
tion for the two main pulsation modes of V391 Peg. The 
349.5 s period is a radial mode: x^(' = 0)=1.5, x^(' — 
1)=53.6, x^(2 < I < 5)>170. The 354.1 s period is a 
dipole mode and again there is only one solution com- 
patible with the data: x^(' = 1)==2.5, x^il = 2)=17.1, 
X^{1 = 0)=47.4, x^(3 < I < 5)>180. These numbers 
translate into a value of the quality-of-fit parameter 
(Press et al. 1986) Q ^ 0.001 for both modes when we 
use an I value different from and 1 respectively. 

4 Discussion 

Thanks to the high quality of the data, this is the first 
time that the mode degree index has been uniquely 
identified from multicolor photometry for the two main 
modes of a star as faint as V391 Peg (V=14.6). To 
our knowledge, conclusive results were obtained only 
for two brighter stars: KPD 2109-1-4401 (V=13.4, Ran- 
dall et al. 2005, see also Jeffery et al. 2004) and Bal- 
loon 090100001, the brightest known sdBV star with 
V=12.1 (Baran et al. 2008, Charpinet et al. 2008). 

The results reported in this article confirm that 
multicolor photometry can set useful identification con- 
straints on the pulsation modes of sdB rapid pulsators, 
provided that the data have a high enough quality. UL- 
TRACAM on a 4 m class (or larger) telescope is an 
ideal instrument for such studies. 
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